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We analyze the separability of the joint state of a collection of two-level systems at finite temper- 
ature T. The fact that only separable states are found in the neighborhood of their thermal equi- 
librium state guarantees that unimpeded thermal decoherence will destroy any initially arranged 
entanglement in a finite time. 
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Significant advances have occurred in the last decade in 
both the mathematical characterization of entanglement 
and in its appHcations, but it remains an important open 
issue how entan^Bment responds to the influence of envi- 
ronmental noise [J . The interaction of a quantum system 
with its surroundings is a generic phenomenon because 
real physical systems cannot be isolated completely from 
external noise. 

So far, it has been remarkably challenging to find any- 
thing generically provable about the process of decoher- 
ence (i.e., disentanglement). For example, if there are 
more than two bodies involved it is still not known how 
to calculate whether a general many-body mixed state 
is separable or not. Here we provide many-body results 
concerning separability of a physically important class of 
mixed states. These arise from consideration of interac- 
tions with thermal reservoirs and properties of thermal 
equilibrium. In many respects these are the least special- 
ized considerations that can be imagined, as every real 
physical system is exposed to a thermal bath at some 
non-zero temperature. 

We suppose only that each sub-system in a many-body 
network of qubits relaxes over time to its thermal equi- 
librium state. We then demonstrate that the many-qubit 
thermal equilibrium state is special, in the sense that it 
is embedded in a finite neighborhood of completely sep- 
arable states. So far as we know, the specification of a 
many-body state with a finite neighborhood of guaran- 
teed separability has been established up to the present 
time only for the identity state (see Zyczkowski, et al. 
[gj). Our result is easily extended to a wide range of 
other physically accessible states, essentially filling all of 
separable state space. This is indicated in Fig. [1] 

Preliminary Considerations:- Our general result is 
established without speculating about approximate or 
"reasonable" measures of many-body entanglement, and 
is first proven abstractly. We then provide new results 
regarding wide-scale entanglement decay in a network of 
two-level systems interacting with independent thermal 
reservoirs. 




FIG. 1: Possible routes of dynamical evolution of states are sug- 
gested by the trajectories sketched in an imaginary non-metrical 
state space where the interior solid line is the boundary between 
separable (S) and non-separable (NS) states. The arrow heads 
simply indicate a direction of evolution. Our first result is that 
evolution to thermal steady state for T > must end at a position 
on the interior of S, as shown for the route labelled T > 0. 



Our results are ultimately based on the fact that the 
thermal equilibrium state itself is separable and diagonal 
in the energy representation. We may represent it by the 
density matrix 



Pq = diag{pi,p2, ...), 



(1) 



with pj > 0, (j — 1, 2, N), for any finite temperature 
T > 0. 

Although our main result is equally valid for an arbi- 
trary number of subsystems, it is convenient in the be- 
ginning to divide the many-body state space Ti. into two 
parts and consider a bipartite system: Ti = H.^ ® HP . 
We assume that both Ti.-^ and TiP are finite many- 
dimensional spaces. An operator p acting on 7i describes 
a state if (1) tr{p) = 1, and (2) tr{pP) > for any pro- 
jection operator P. That is, p is a positive Hermitian 
operator. 

We will work from the fact that a bipartite density ma- 
trix p is non-separable (entangled) if and only if there ex- 
ists a Hermitian "entanglement witness" matrix W with 
the two properties Q: 

tr{Wp) < 0, and (2) 
tr{W(7) > 0, for all separable states a. (3) 
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Finite Neighborhood of Separability:- We assert that 
for a density matrix with the thermal equilibrium form 
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([T]) there exists a neighborhood U of po such that all 
states in that neighborhood are separable. The result of 
Zyczkowski, et al. can be regarded as the limiting 
case of our result as T — cx3. We will prove our assertion 
by showing that the opposite case must be false - i.e., 
the assumption that an entangled state exists in every 
neighborhood of po leads to a contradiction. 

First suppose that a neighborhood C/„ of po has an en- 
tangled state Pn in it. Then consider an infinite sequence 
of smaller neighborhoods, C/„+i C C/„, all having at least 
one entangled state. This implies a sequence of entangled 
states Pn converging to po as n ^ oo. 

According to ([2]) and ([3]), for each entangled state p„ 
there must exist a witness operator W„ that satisfies both 



triWnPn) < 0, and tr{Wn(T) > 0, 



(4) 



for Pn and all separable states a. Therefore, for the 
sequence of states and witnesses, we get (as n oo): 
tr{WnPn) tr{Wpo),tr(Wna) tr{Wa), for all sepa- 
rable states cr. Moreover, by combining these expressions, 
and remembering the separable character of po, we have 



tr{Wpo) < 0, 
tr{Wa) > 0. 



(5) 
(6) 



Then since po is acceptable as a possible a, we conclude 
that satisfying both ^ and ([6]) requires 



tr{Wpo) ^piWii +P2W22 



= 0. 



(7) 



At this point a useful observation is that the state rep- 
resented by (Ti = c?ia(7(l, 0, 0)yi c?m(7(l, 0, 0)b is 
obviously separable, and if W is an entanglement witness, 
then tr{Wai) > requires Wn > 0, and by extension 
Wjj > for all j. Thus, since we are concentrating on 
the case in which the pjS in ^ are all greater than zero, 
we conclude from ([7]) that all the diagonal elements Wjj 
must vanish. Therefore, we have trW — tr{WI) — 0, 
where / is the identity matrix. Since W is not a zero 
matrix, there exists a product projection P (E) Q such 
that tr{W P iSi Q) 7^ (see, e.g., Q). By adding and sub- 
tracting / and using tr{WI) = again, this is sufficient 
to show that 



tr[W{P (g) Q)] = -tr[W{I - P ® Q)] 7^ 0. 



(8) 



Both P ®Q and [I ~ P ® Q) are (at least proportional 
to) separable states, so ^ cannot hold for both. This is 
the needed contradiction that proves our main result. 

An immediate implication of the existence of a sepa- 
rable neighborhood for a system consisting of M qubits, 
each coupled to a local thermal heat bath, is that unim- 
peded thermal decoherence must destroy the entangle- 
ment of every initial state in a finite time, an effect we dis- 
cussed for zero-temperature two-atom spontaneous emis- 
sion in a previous note (Q, referred to below as YE for 
short). This has been labelled ESD for early-stage deco- 
herence or "entanglement sudden death" , indicating that 



in order to reach po asymptotically the system must en- 
ter the separable neighborhood of po after only a finite 
time. 

M-Qubit Systems Under Thermal Noise:- The dynam- 
ics of pure thermal decoherence is completely determined 
by the reduced density matrix of the M-qubit system, ob- 
tained by tracing over the other qubits and the thermal 
reservoir's variables. In the familiar Born-Markov ap- 
proximation, when each qubit is in contact with a broad- 
band harmonic reservoir at temperature T, one finds a 
compact Lindblad master equation (e.g., see Q) for M 
qubits {h = 1): 



— p = -i[Hsys,p]+ C{p), 



(9) 



where H. 
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superoperator is C{p) = 
\pC,A, - \G.fi%P) and 
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and the Liouvillian 

the Lindblad operators 
Ci2 = V^a+, = 



n -\- l)rtT_, Cii — \/^T(T_|_. Here n is the mean num- 
ber of thermal reservoir quanta: n = Yjifi^l^^'^ — !)■ 

The general solutions of the master equation © can 
be expressed compactly by the Kraus operators Kf 
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operators Kj 



(10) 



with 



elementary Kraus 
j^^^^j^^j — I are the same for each qubit, given 
equal temperatures for all the reservoirs, and are explic- 
itly given by 
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(11) 

(12) 

(13) 
(14) 



where the time-dependent Kraus matrix elements are 
7(t) = exp [-^T{2n + l)t] , w(t) = ^l-j'^{t). The op- 
erators Ki , K2 provide the transitions from the excited 
state |-|-) to the ground state |— ) caused by stimulated 
and spontaneous emission whereas the operators K^, K4 
account for absorptive transitions from the ground state 
| — ) to the excited state |-|-). 

For the Af-qubit system described by the master equa- 
tion above, the steady state will be given by ([T]) . Let 
p{t) denote the state of M qubits at t. Then p{t) po 
as t 00. Our generic result implies that a finite ESD 
time, ^csd; exists for an arbitrary initial state. A simple 
expression for ^osd is given in YE for two-atom sponta- 
neous emission. 
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Pair Entanglement in Thermal States:- For a two- 
qubit subsystem, our generic results established above 
can be realized in a more concrete way by explicitly 
solving the time evolution. Bipartite dynamics under 
thermal noise for T > has been treated previously in 
considering various aspects of entanglement, e.g., gener- 
ation^ I, fragility [sl, ^ , influence of squeezed reservoirs 
p^.[lir 1^, Brownian particle diffusion 13], and univer- 
sality of BSD [3, [l^ . We now show explicitly that when 
each qubit relaxes to its equilibrium state asymptotically, 
the entanglement between any pair of qubits will always 
terminate in a finite time irrespective of the initial states 
of the qubits. 

For any qubit pair, the steady state will be the diagonal 
matrix ^ with iV = 4: 



Pst = dmg(pi,p2,P3,P4), 



(15) 



in the standard basis | + +), | H — ), | — h), | ), where 

the thermal probabilities are pi = n^/(2n -f 1)^, P2 — 
P3 = n{n + l)/{2n + 1)^, pi = {n+ 1) V(2n + 1)^. 

A standard measure of entanglement for our two-qubit 
system is Wootters' concurrence [l^, denoted C. By 
construction, the concurrence varies from C — for a 
separable state to C = 1 for a maximally entangled state. 
For the two qubits A and B we have: 



C^-^ (p) = max [0, A(p)] , where 
A(p) = - - \/h - 



(16) 
(17) 

and the Xi are the eigenvalues in decreasing order of the 
matrix = p(cr^ (g) ay )p*{ay ® Uy). Here p* denotes 
the complex conjugation of p in the standard basis given 
above and cr^'^ are the Pauli matrices expressed in the 
same basis. 

Formulated in this way, a necessary and sufficient con- 
dition for p to be separable (zero entanglement) is A < 0. 
We have noted previously ^\ that the negative values 
of A possess quantum state information not contained 
in the concurrence C{p). Both A = and A < sig- 
nal that the two-qubit density matrix p is separable, and 
when A(p) < 0, we call the state p a super-separable 
state. The distinction for quantum state trajectories is 
shown in Fig. [Jl where some trajectories terminate in- 
side the S zone (A < 0) and some terminate exactly 
at the edge of the S zone (A — 0). For the ther- 
mal steady state pst, straightforward calculations yield 
A(p,t) = -2VP2P3 = -2n{n + l)/(2n + \f < 0, so the 
steady state is separable for any reservoir temperature 
including zero, but super-separable for all T > 0. 

In order for the system to reach the zone of super- 
separable states, A must go from positive to negative 
values. Since A is a real number, this means it must cross 
the value A = 0, and this suggests a different graphical 
representation of trajectories with metrical elements that 
can assist interpretation, as in Fig. [2l Whenever the 
long-time steady state of a two-qubit system is super- 
separable, the system must experience ESD. This means 
that it must become separable after only a finite time. 



as in the lower two curves in the figure, where the dots 
identify two locations for t^sd- This possibility was ap- 
parently first noticed in several different physical contexts 
independently (see 0, [H [H, [H, [H, Ho] ) . A two-qubit 
ESD time was evaluated in YE, and ESD was first re- 
ported experimentally by the Davidovich group [IH and 
also recently by the Kimble group in atomic ensembles 

In some cases of interest the ESD times t^sd can be 
determined explicitly. It is instructive to follow en- 
tanglement evolution for the two Bell states |^'±) = 
^(| + -)±|-+>). From (Uni), the evolving BeU density 
matrix at t is given by 



P{t) = 



a{t) 

b{t) z{t) 

z{t) c{t) 

d{t) 



(18) 



where the time-dependent matrix elements are given by 
the following: a(t) = p2'^^ +P37^i^^, b{t) = c(t) = pij'^ + 
P2l^+P3{l^ + l+UJ^),d{t) ^piuj^+p3j^u;^,z{t) ^±\l^. 
Then it is straightforward to compute A(p(t)) = 2|z(i)| — 
1^a(t)d{t). Thus, we conclude that when 



i > tcsd = Y 



l + 2^plp? 



(19) 



the initial Bell states |^'±) become completely disentan- 
gled. As suggested by the two solid curves in Fig. [21 
the higher the temperature the shorter the disentangle- 
ment decay time t^sd ■ The upper solid curve is consistent 
with YE, showing no Bell disentanglement in a finite time 
when the environments are at zero temperature. 

Conclusion:- The central element of this note is our 
demonstration of a finite neighborhood of separable 
states surrounding the thermal equilibrium state, but the 
proof could equally well have been given for any state in 
a finite dimensional space having diagonal form and pos- 
itive definite elements. In fact continuously many other 



A(co) 




CO(t=oo) 



FIG. 2: Four types of disentanglement trajectory are plotted via 
their A values, where the horizontal axis covers t = [0, oo], scaled 
according to Lo{t) given below Eqn. ||14|I . The solid lines show 
evolution for T = and T > and each is appropriate to both Bell 
states, which start with A = 1. While ESD is universal for T > 0, 
the dashed lines, starting from A = 1/2, show evolution of mixed 
states and demonstrate that ESD can also occur for T = 0. The 
dots where two curves cross the lj axis identify tegj times. 
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states of this type are also physically realizable, in some 
cases easily. Any completely incoherent partially excited 
state would suit, and an obvious example is the state 
associated in spin resonance with finite but "negative" 
temperature or in laser physics with a positive partial 
inversion. 

All of these also have finite neighborhoods of separabil- 
ity and all such neighborhoods are clearly linked to each 
other. It is an open and interesting question whether the 
volume of this set of states can make up a significant or 
even finite fraction of the total separable volume. We 
note that our demonstration did not require knowing the 
finite number M of qubits under consideration, and so 
the existence of a finite separable neighborhood is inde- 
pendent of M . The result is topological, and does not 
confiict either with the remark of Zyczkowski, et at, 
that the volume of separable states decreases with system 
size, or the proof by Eisert, et at, [23!| that non-separable 
states are dense in continuously infinite Hilbert spaces. 

Not surprisingly, ESD and other forms of dissipative 
entanglement dynamics and their control are of current 



interest, with examinations reported of many different 
dynamical processes and realizations of qubits, includ- 
ing coupled mechanical oscillators llSll . cubits in a spin 
chain and coupled to an Ising chain 'M, 25] , multi-cavity 
QED [2E[27[, spin ensembles coupled via lossy photonic 
channels 23 j qubit-qutrit combination I29l. two-qubit 
decoherence dynamics [13, llll, HI, HI, l34j and multiple 
noises [35], to name a few. 

To summarize, the main purpose of this note is to ex- 
tend the study of entanglement in thermal and thermal- 
like environments to general multipartite systems with 
arbitrary initial states. We have achieved this goal by 
exploring topological properties of the state space of a 
many-body system. As an illustration, we have demon- 
strated the necessary thermal decoherence of two-qubit 
systems. The suppression of thermal decoherence and 
treatment of interacting qubits will be considered in fu- 
ture publications. 
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